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Abstract 

It is shown that for the N-neighbor and K-state cellular automata, the class II, 
class III and class IV patterns coexist at least in the range < A < 1 — The 
mechanism which determines the difference between the pattern classes at a fixed A 
is found, and it is studied quantitatively by introducing a new parameter F. Using 
the parameter F and A, the phase diagram of cellular automata is obtained for 
5-neighbor and 4-state cellular automata. 

PACS: 89.75.-k Complex Systems 



1 Introduction 



Cellular automata (CA) has been one of the most studied fields in the re- 
search of complex systems. Various patterns has been generated by choosing 
the rule tables. Wolfram [1] has classified these patterns into four rough cat- 
egories: class I (homogeneous), class II (periodic), class III (chaos) and class 
IV (edge of chaos). The class IV patterns have been the most interesting 
target for the study of CA, because it provides us with an example of the 
self-organization in a simple system and it is argued that the possibility of 
computation is realized by the complexity at the edge of chaos[2,3,4]. 

A more detailed classification of CA, has been carried out mainly for the 
elementary cellular automata (3-neighbor and 2-state CA)[5,6], in which the 
pattern is studied quite accurately for each rule table. However the number 
of rule table in N-neighbor and K-state cellular automata grows like K K . 
Therefore except for a few smallest combinations of the N and K, the num- 
bers of the rule tables become so large that studies of the CA dynamics for 
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all rule tables are impossible even with the fastest supercomputers. 

On the other hand the rule table of edge of chaos is rare in the whole CA 
rule table space, therefore it is important to find a set of parameters by which 
the pattern classes could be classified, and to determine a phase diagram of 
CA in these parameter space, even if it is a qualitative one. 

Langton has introduced A parameter and argued that as A increases the 
pattern class changes from class I to class II and then to class III. And 
in many cases, class IV behavior is observed between class II and class III 
patterns[7,3,8]. The A parameter represents rough behavior of CA in the rule 
table space, but finally does not sufficiently classify the quantitative behavior 
of CA. It is well known that different pattern classes coexist at the same A. 
Which of these pattern classes is chosen, depends on the random number. The 
reason or mechanism for this is not yet known; we have no way to control the 
pattern classes at fixed A. And the transitions between a periodic to chaotic 
pattern classes are observed in a rather wide range of A. In Ref.[8], a schematic 
phase-diagram was sketched. However a vertical axis was not specified. There- 
fore, it is has been thought that more parameters are necessary to arrive at a 
more quantitative understanding of the rule table space of the CA. 

In this article, we will report a mechanism which discriminates the pattern 
classes at a fixed A. The mechanism is closely related to the structure of the 
rule tables and is expressed by the numbers of rules which breaks strings of 
quiescent state. For N-neighbor and K-state CA, it is found that in the region 
1/K < X < 1 — 1/K, the class II, class III and class IV pattern classes coexist. 

This property is studied quantitatively by introducing a new parameter F, 
which is taken to be orthogonal to A. In the region 1/K < X < 1 — 1/K, the 
maximum of F correspond to class III rule tables while minimum of F, to 
class II or class I rule tables. Therefore the transition of the pattern classes 
takes place somewhere between these two limits without fail. If we determine 
the region of F, where the transitions of the pattern classes take place, we 
could obtain the phase diagram in A — F plane. 

The determination of the phase diagram is carried out for 5-neighbor and 
4-state CA. In this case, phase boundary is not sharp but has some range in 
A — F plane. The region has a gentle slope as a function of A, and extends 
over the range 0.2 < A < 0.8. This means that for this CA the edge of chaos 
could be found at least in this range in A. 

In section 2, we will briefly summarize our notations and present a key dis- 
covery, which leads us to the understanding of the structure of the rule table 
and pattern classes. It strongly suggested that the rules which break strings 
of the quiescent states play an important role for the pattern classes. 

In section 3, we classify rule tables according to the destruction and con- 
struction of strings of the quiescent states, and carry out the replacements 
of the rules to change the chaotic pattern class into periodic one and vice 
versa while keeping A fixed. The reason why the patten classes changes by the 
replacements is discussed, and we will show that in the region 1/K < X < 
1 — 1/K, the change of the pattern classes takes place without fail by the 
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replacements. 

In section 4, the result obtained in section 3 is studied quantitatively by 
introducing a new parameter F. Using F and A, we determine the phase dia- 
gram in the A — F plane for 5-neighbor and 4-state CA. 

Section 5 is devoted to conclusions and discussions. 



2 Summary of CA and a key discovery 



2.1 Summary of Cellular Automata 



In order to make our arguments concrete, we focus mainly on the one- 
dimensional 5-neighbor and 4-state CA in the following, however, the qual- 
itative conclusions hold true for other CAs. This point will be discussed in 
subsections 3.2. 

We will briefly summarize our notation of CA[1,3]. In our study, the site 
consists of 150 cells having the periodic boundary condition. The states are 
denoted as s(t, i). The t represents the time step which takes an integer value, 
and the % is the position of cells which range from to 149. The s(t,i) takes 
values 0, 1, 2, and 3, and the state is taken to be the quiescent state. The 
set of the states s(t,i) at the same t is called the configuration. 

The configuration at time t + 1 is determined by that of time t by using 
following local relation, 

s(t + l,i)= T(s(t, i-2), s(t, s(t, i), s(t, i + 1), s(t, i + 2)). (1) 



The set of the mappings 

T(fi, v, k, p, o) = V, (», etc. = 0, 1, 2, 3) (2) 



is called the rule table. The rule table consists of 4 5 mappings, which are 
selected from a total of 4 1024 elements. 
The A parameter is defined as [3] 



where iVj, is the number in which 77 in Eq.2 is not equal to 0. In other words 
the A is the probability that the rules do not select the quiescent state in next 
time step. In the following we set the rule tables randomly with the probability 
A. We choose 1024 — N h rules randomly, and set r\ = in the right hand side of 
Eq. 3. For the rest of the N h rules, the 77 picks up the number 1, 2, 3 randomly. 
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The initial configurations are also set randomly. 

The time sequence of the configurations is called a pattern. The patterns 
are classified roughly into four classes established by Wolframfl]. It has been 
known that as the A increases the most frequently generated patterns change 
from homogeneous (class I) to periodic (class II) and then to chaotic (class 

III) , and at the region between class II and class III, the edge of chaos (class 

IV) is located. 



2.2 Correlation between pattern classes and rules which break strings of qui- 
escent state at a fixed A 

In order to find the reason why the different pattern classes are generated 
with the same A, we have started to collect rule tables of different pattern 
classes, and tried to find the differences between them. We have fixed at A = 
0.44 (Nh = 450), because at this point the chaotic, edge of chaos, and periodic 
patterns are generated with a similar ratio. By changing the random number, 
we have gathered a few tens of the rule tables and classified them into chaotic, 
edge of chaos, and periodic ones. 

In this article, a pattern is considered the edge of chaos when its transient 
length [3] is longer than 3000 time steps. 

First, we study whether or not the pattern classes are sensitive to the initial 
configurations. We fix the rule table and change the initial configurations. The 
details of the patterns depend on the initial configurations, but the pattern 
classes are not changedfl]. Thus the difference of the pattern classes is due to 
the differences in the rule tables, and the target of our inquiry has to do with 
the differences between them. 

For a little while, we do not impose a quiescent condition (QC)[3], 
T(0, 0, 0, 0, 0) = 0, because without this condition, the structure of the rule 
table becomes more transparent. This point will be discussed at the footnote 
4 in section 4. 

After some trial and error, we have found a strong correlation between the 
pattern classes and the QC. For class II patterns, the probability of the rule 
table, which satisfies the QC is much larger than that of the class III patterns. 
This correlation has suggested that the rule T(0, 0, 0, 0,0) — h, h ^ 0, which 
breaks the string of the quiescent states with length 5, pushes the pattern 
toward chaos. We anticipate that the similar situation will hold for the strings 
of quiescent states with length 4. 

We go back to the usual definitions of CA. In the following we discuss CA 
under QC, T(0, 0, 0, 0, 0) = 0. We study the correlation between the number 
of the rules of Eq.4 and the pattern classes: 

T(0,0,0,0,i) = h, T(i, 0, 0, 0, 0) = h, (i, h = 1, 2, 3). (4) 
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Fig. 1. The pattern classes at A = 0.44. The quiescent state is shown by white dot, 
while other states are indicated by black point. Fig. 1(a) corresponds to N 4 = 4, 
Fig.l(b), Fig.l(c), and Fig.l(d), to N 4 = 3, Fig.l(e), Fig.l(f), and Fig.l(g), to 
N4 = 2, and Fig. 1(h) corresponds to N4 = 1. 



These rules break length 4 strings of the quiescent states, and will also push 
the pattern toward chaos We denote the total number of rules of Eq.4 in 
a rule table as N4. We have collected 30 rule tables and grouped them by the 
number N4. We have 4 rule tables with N± > 4 , 13 rule tables with N4 = 3, 
9 rule tables with N 4 = 2 and 4 rule tables with iV 4 < 1. When N 4 > 4, all 
rule tables generate chaotic patterns, while when N 4 < 1, only periodic ones 
are generated. At N4 = 3 and N 4 = 2, chaotic, edge of chaos, and periodic 
patterns coexist. Examples are shown in the Fig.l. The coexistence of three 
pattern classes at N4 = 3 is seen in Fig. 1(b), Fig. 1(c) and Fig. 1(d) and that 
of N4 = 2 is exhibited in Fig. 1(e), Fig. 1(f) and Fig. 1(g). 

The strong correlation between N4 and the pattern classes has been ob- 
served in this case too, as anticipated. These discoveries have provided us with 
a key hint leading us to the hypothesis that the rules, which break strings of 
the quiescent states, will play a major role in the pattern classes. 



2 Similar ideas had been noticed by Wolfram and Suzudo with the arguments of 
the unbounded growth[l] and expandability [9]. 
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3 Structure of rule table and pattern classes 

3.1 Structure of rule table and replacement experiment 

In order to test the hypothesis of the previous section, we classify the rules 
into four groups according to the operation on strings of the quiescent states. 
In the following, Greek characters in the rules represent groups 0, 1, 2, 3 while 
Roman, represent groups 1, 2, 3. 
Group 1: T(p, z/, 0, p, a) = h. 

The rules in this group break strings of the quiescent states. 

Group 2: T(/i, z/, 0, p, a) = 0. 

The rules of this group conserve them. 

Group 3: T(fi, z/, i, p, a) = 0. 

The rules of this group develop them. 

Group 4: T(/z, z/, i, p, a) = I. 

The rules in this group do not affect string of quiescent states in next time 
step. 

The sum of the numbers of the group 1 and group 2 rules is 256, while 
that of group 3 and group 4 rules is 768. The number of each group of rules 
included in the rule table is determined mainly by the probability A, therefore 
it suffers from fluctuation due to randomness. 

The group 1 rules are further classified into five types according to the length 
of string of quiescent states, which they break. These are shown in Table 1. 

Table 1 

The classification of the rules in group 1. 



type 


Total Number 


Name 


Replacement 


T(0,0,0,0,0) = h 


1 


D5 


RP5,RC5 


T(0,0,0,0,i) = h 


3 


D4 


RP4,RC4 


T(i,0, 0,0,0) = h 


3 






r(0,0,0,i,a) = h 


12 






0, 0, 0, m) = h 


9 


D3 


RP3,RC3 


T(n,j,0,0,0)=h 


12 






T{n, j,0,0,m) = h 


36 


D2 


RP2,RC2 


T(i, 0,0, 1, a) = h 


36 






T(n,j,0,l,a) = h 


144 


Dl 


RP1,RC1 



The D5 rule is always excluded from rule tables by the quiescent condition. 

Our hypothesis presented at the end of the section 2 is expressed more 
quantitatively as follows; the numbers of the D4, D3, D2, and Dl rules shown 
in Table 1 will mainly determine the pattern classes. 

In order to test this hypothesis we artificially change the numbers of these 
rules in Table 1 while keeping the A fixed. For D4 rules, we carry out the 
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replacements defined by the following equations, 

T(0, 0, 0, 0, i) = h -> T(0, 0, 0, 0, i) = 0, 
or T(i, 0, 0, 0, 0) = h -> T(i, 0, 0, 0, 0) = 0, 
T(n, u, j, p, a) = -> T(p, i/, j, p, a) = /, 



where except for h, the groups p, z/, p, a, j and / are selected randomly. 
Similarly the replacements are generalized for D3, D2, and Dl rules , which 
are denoted as RP4 to RP1 in Table 1. They change the rules of group 1 to 
that of group 2 together with group 3 to group 4 and are expected to push 
the rule table toward the periodic direction. 
The reverse replacements for D4 are 

T(0, 0, 0, 0, i) = -> T(0, 0, 0, 0, i) = h, 

or T(i, 0, 0, 0, 0) = -> T(i, 0, 0, 0, 0) = h, , . 

^(P, ^, J, p, <r) = Z -> T(p, z/, j, p, <r) = 0, ^ ' 



which will push the rule table toward the chaotic direction. In this case, the 
groups h, p, v, j, p, and a are selected randomly. Similarly we introduce the 
replacements for D3, D2, and Dl, which will be called RC4 to RC1 in the 
following. 

By the replacement of RP4 to RP1 or RC4 to RC1, we change the numbers 
of tfhe rules in Table 1 while keeping the A fixed. We denote these numbers 
A^ 4 , A^ 3 , A^ 2 , and N\ for D4, D3, D2, and Dl rules, respectively. The examples 
of the replacement experiments are shown in Fig.2. 

The rule table of Fig. 2(a) is obtained randomly with probability A = 0.6. 

Table 2 

The numbers of the rules and the replacements for each figures shown in Fig.2 



Figure 


7V 4 


N 3 


N 2 




RP4 


RP3 


RP2 


RP1 


Fig.2(a) 


3 


22 


53 


96 














Fig.2(b) 





14 


53 


96 


3 


8 








Fig.2(c) 





13 


53 


96 


3 


9 








Fig.2(d) 





12 


53 


96 


3 


10 








Fig.2(e) 


1 


7 


53 


96 


2 


15 








Fig.2(f) 


1 


6 


53 


96 


2 


16 








Fig.2(g) 


2 


1 


53 


96 


1 


21 








Fig.2(h) 


2 





53 


96 


1 


22 









At A = 0.6 most of the randomly created rule tables generate chaotic patterns. 
Fig.2 (b) is obtained by the 3 RP4s and 8 RP3s, and the numbers of the rules 
become A^ 4 = and iV 3 = 14, respectively. These numbers are summarized 
in Table 2. By this replacements, Fig.2 (b) shows an edge of chaos behavior. 
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Fig. 2. Example of the replacement experiments at A = 0.6. Fig. 2 (a) is obtained 
randomly with probability A = 0.6, by the method explained in subsection 2.1. 
Fig.2(b) to Fig.2(h) are obtained by the replacements of the rules of Fig.2(a), which 
are summarized in Table 2. 



When N 4 = 0, until N 3 = 22 to N 3 = 15, the rule tables generate chaotic 
patterns. Fig. 2(c) is obtained by one more RP3 replacements for Fig. 2(b) rule 
table. It shows a periodic pattern with a rather long transient length. One 
more replacement of RP3 for the Fig. 2(c) rule table is shown in Fig. 2(d), 
where the transient length becomes shorter. 

Similar replacement experiments for JV4 = 1 and N4 = 2 cases are shown 
in Fig. 2(e), Fig.2(f), and Fig. 2(g), Fig. 2(h), respectively. In these cases, until 
N 3 > 8 and N 3 > 2, rule tables generate class III patterns at A^ 4 = 1 and 
N 2 = 2, respectively. In these examples, edges of chaos are observed between 
classes III and II behaviors. 

For each A = 0.8, 0.75, 0.7, 0.6, 0.5, 0.4, 0.3 and 0.2 point, we have carried 
out a several hundreds to a few thousands replacement experiments. In these 
replacements, we have succeeded in changing chaotic rule tables to a periodic 
ones by the replacements of RP4 to RP2, and vice versa by RC4 to RC2. For 
the replacements RP4, RP3 and RP2, we have observed no example that the 
rule table moves to chaotic direction. The converse is true for RC4, RC3 and 
RC2. Exceptions are observed only in the replacements RP1 and RC1, which 
will be discussed in section 4. 
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3.2 Chaotic and periodic limit at fixed A in N-neighbor and K-state CA 

Let us study the effects of the replacements theoretically in the N-neighbor 
and K-state CA. In this general case too, the rule tables are classified into four 
groups. We denote the number of the group 1 rules as N(gl), and similarly 
for the numbers of the other groups. These numbers satisfy the following sum 
rules. 

N{g2) + N(g3) = K N (1 - A), N(gl) + N(g4) = K N X, 

N(gl) + N(g2) = K N ~ l , N(g3) + N(gA) = K N ~ l (K — 1). {) 

The individual numbers N(gi) are determined by the probability A. They 
are summarized in the Table 3. They suffer from fluctuations due to random 
number, however in this subsection, we neglect the fluctuations. 

The replacements to decrease the number of the group 1 rule while keeping 

Table 3 

The classification of N-neighbor and K-state CA rules into four groups. The fii 
represent to K-l and h,i and /, 1 to K-l. 





rule 


N(gi) 


group 1 


T(ui,u 2 , ...,0, ...,u N ) = h 




group 2 


r(/ii,^ 2 ,-,o,...,/ijv) = o 


K"-'(l-\) 


group 3 


T(ui,u 2 , ...,u N ) = 


K N ~\K -1)(1-X) 


group 4 


T(m,n2, ...,/xjv) = i 


K N ~ L (K — 1)A 



the A fixed are given by, 

N(gl) -> N(gl) - 1, N(g2) -> N(g2) + 1 

N(g3) N(g3) - 1, iV(^4) ^ 7V(y4) + 1 (8) 

In the case of 5-neighbor and 4-state CA, they correspond to RP4 to RP1. 

These replacements stop either when N(gl) = or N(g3) = is reached. 
Therefore when N(gl) < N(g3), which corresponds to A < (1 — j^) in A, all 
the group 1 rules are replaced by the group 2 rules. In this limit, quiescent 
states at time t will never be changed, because there is no rule which converts 
them to other states, while the group 3 rules have a chance to create a new 
quiescent state in the next time step. Thus the number of quiescent states at 
time t is a non- decreasing function of t; therefore, the pattern class should 
be class I (homogeneous) or class II (periodic), which we call periodic limit. 
Therefore the replacements of Eq.8 push the rule table toward the periodic 
limit. 

Let us discuss the reverse replacements of Eq.8. In these replacements, if 
N(g2) < N(gA), all group 2 rules are replaced by the group 1 rules, except for 
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the quiescent condition. In this extreme reverse case, all the quiescent states 
at time t are converted to other states in next time step, while group 3 rules 
create them at different places. Then this will most probably develop into 
chaotic patterns. This limit will be called chaotic limit, which is reached in 
the region ■h < A. We should like to say that atypical rule table and initial 
condition might generate a periodic patterns even in this limit. But in this 
article, these exceptional cases are neglected. 
Therefore in the following region, 

i<A<l-l (9) 



all the rule tables are located between these two limit, and by the replacements 
of Eq. 8 and their reverse ones, the changes of the pattern classes take place 
without fail. This explains the validity of the hypothesis of previous section. 
And we have found a method to control the pattern classes at fixed A. 



4 Phase diagram of 5-neighbor and 4-state CA in A-F plane 



In the previous section, we have found that rule table is located somewhere 
between chaotic limit and periodic limit, in the region ^ < A < 1 — i. In 
order to express the position of the rule table quantitatively, we introduce 
new parameter F, which provides us with a new axis (F-axis) orthogonal to 
A. Minimum of F is the periodic limit, while maximum of it corresponds to 
chaotic limit. In this section, we determine the parameter F, for 5-neighbor 
and 4-state CA. 

As a first approximation, the parameter F is taken to be be a function of 
the numbers of the rules D4, D3, D2 and Dl, which have been denoted as A4, 
A 3 , A 2 and Ni, respectively. We proceed to determine F(A 4 , A 3 , A 2 , Ni) by 
applying simplest approximations and assumptions 

We have observed in the replacement experiments, that the position of the 
rule table in F moves toward chaotic direction, when A4 or A3 or N2 increases. 
Examples are shown in Fig. 2 and Table 2. However for Dl, replacements RP1 
and RC1 sometimes look like random walk on F-axis, around the region where 
the transition of the pattern class is taken placef 3 } Therefore F will be a com- 
plicated function of Ni and determination of it will be difficult. 

However the number jVj is rather large, therefore we apply mean field ap- 
proximation for Ai. We replace Ai by its average, and measure F from this 

3 The effect of Dl rule is to change an isolates quiescent state to other states in the 
next time step. This effects may easily be compensated by the creation of quiescent 
states by group 3 rules in one time step. This may be a reason that the replacements 
of RP1 and RC1 some times look like random walk. 
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background; namely the A"i dependence is dropped from F and take F = at 
N 4 = N 3 = N 2 = 0. 

Here, we apply Taylor series expansion for F at this point, and approximate 
it by the linear terms in A4, iV 3 and N 2 . 



F(N 4 , N 3 , N 2 ) ~ c 4 N 4 + c 3 N 3 + c 2 N 2 . (10) 

where c 4 = dF/dN^, similar for c 3 and c 2 . They represent the strength of 
the effects of the rules D4, D3, and D2 to push the rule table toward chaotic 
direction. This definition is symbolic, because N4 is discrete. 

The measure in the F is still arbitrary. We fix it in the unit where the 
increase in one unit of N4 results in the change of F in one unit. This corre- 
sponds to divide F in Eq.10 by c 4 , and to express it by the ratio c 3 /c 4 (r 3 ) 
and c 2 /c 4 (r 2 ). 

Before we proceed to determine r 3 and r 2 , let us interpret the parameter F 
geometrically. Most generally, the rule tables are classified in 1024-dimensional 
space in this CA. The rule tables, at the boundary of the class III and class 
II pattern classes at fixed A, form a hyper-surface in this space. We map the 
points on hyper-surface into 3-dimensional (jV 4 , N 3 , N 2 ) space. They will be 
located in some region in the 3-dimensional space. We introduce a surface 
F(A" 4 , A3, A2) = $ in order to line up these points. F-axis is a normal line of 
the surface F(A" 4 , A3, N 2 ) = $. In Eq.10, we approximate it by a plane. 

Our strategy to determine r 3 and r 2 is to find the regression plane in JV4, 
N 3 , and A" 2 space. It is equivalent to fix the F-axis in such a way that the 
projection of the distribution of transition points on F-axis, (F crit ) looks as 
narrow as possible. The quality of our approximations and assumptions re- 
flects the width of the distribution of F cr i t . 

In the least square method, our problem is formulated to find r 3 and r 2 , 
which minimize the quantity, 

s(r 3 ,r 2 ) = \ Y,( F UNlNiNi) - FUNi,Ni,Ni)) 2 , (11) 

where % and j label the rule tables on the hyper-surface mapped into 3- 
dimensional space. We solve the equations, dS/dr 3 = and dS/dr 2 = 0, 
which are 

rs £ W) 2 + r 2 £ SNpSN? = - £ SNpSN*, 

r 3 £ SNpSN? + r 2 £ W) 2 £ SNpSN?, (12) 

i,j i,j i,j 
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where 5iVj j = N\ - N[, similar for 5N l 3 ' j and 5N l 2 ' j . 

By artificially carrying out the replacements of RP4, RP3 and RP2 or RC4, 
RC3 and RC2, we look for critical combinations of N±, N 3 , and N 2 , at which 
the change of the pattern classes are observed. The examples of the critical 
combinations are presented in the lines of Fig. 2(b), Fig. 2(e), and Fig.2(g) of 
Table 2. 

The numbers of the critical combinations (Nc tot ), which are used to de- 
termine r 3 and r 2 in Eq. 12, are summarized in the Table 4. The r 3 and r 2 
are determined for each A. Their results are also shown in Table 4, where the 
errors are estimated by the jackknife method. 

The results show that the coefficients are positive, and satisfy the order, 

c 4 > c 3 > c 2 . (13) 

It means that the effects to move the rule table toward chaotic limit on the 
F-axis are stronger for the rules which break longer strings of the quiescent 
states. 

The order in Eq.13 is understood by the following intuitive arguments. If six 
D4 rules are included in the rule table, the string of the quiescent states with 
length 5 will not develop. Similarly, if 33 D3 rules are present in the rule table, 
no length 4 string of the quiescent states could be made. These are roughly 
similar situations for pattern classes. Thus the strength of the D3 rules will 
be roughly equal to 6/33 of that of D4 rules, similar for the strength of the 
D2 and Dl rules Q. 

Using these results for r 3 and r 2 , we calculate Fs for each critical combi- 
nation (F cr it). For A = 0.6, they are shown in Fig.3a for 101 data points. The 
number distributions of them are displayed in Fig. 3b. From these 101 data, 
we calculate average, < F crit > and standard deviation o of the distribution 
of F C ri t . Similar calculations are carried out for other A points. The results 
for < F crit >, < Fcrit > ±cr and maximum of F (F max ) determined by the 
coefficients r 3 and r 2 are shown in the Fig. 4. This is a phase diagram in A — F 
plane. 

The phase diagram is natural in the sense that as A increases, the propor- 
tion of the chaotic rule table region increases in the total F range. The two 
lines < Fcrit > ±c in the figure, indicate that about 68% of the critical com- 
binations are located in this region, if the normal distribution is assumed for 

4 If we do not impose the quiescent condition, Eq.13 becomes 

c 5 > c 4 > c 3 > c 2 . 

Therefore the correlation between pattern classes and the existence of D5 rule is 
stronger than the correlation between those and the number of D4 rules. If we start 
our study within the quiescent condition, we may make a longer detour to find the 
hypothesis of section 2 and get the qualitative conclusion of section 3. 
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Table 4 

The number of the critical combinations Nc tot and the coefficients and ri. 



A 




Relative Strength 


Error 


Nc tot 


n 9 


' 3 


U . 1UOO 


n no90 


± ± ij 




' z 


0.0153 


0.0025 




u.o 


r 3 


U. 1UOO 
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' z 
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0.0019 
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0.1631 


0.0063 
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0.0042 




0.75 


^3 


0.1505 


0.0114 
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0.0419 


0.0065 
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0.1315 


0.0168 
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T2 
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0.0038 
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Fig. 3. Distribution of critical combinations F cr u at A = 0.6. Fig. 3a shows 101 indi- 
vidual F cr a points. We label these 101 points by Nc. In Fig. 3b, number distribution 
of the critical combinations in the bin is displayed. The region of distribution of 
F cr it is divided into 6 bins. 

Fcrit- We call this region as transition region. We think that the class II and 
class III patterns would not be separated by a line in A — F plane. Because 
the transition region in Fig. 4 is a mapping of the hyper-surface, which sepa- 
rates the two pattern classes at a fixed A in 1024 dimensional space, into the 
normal line of the regression plane in the 3-dimensional space. Therefore the 
two parameters will not be enough to draw the boundary by a line. 

However we should like to notice that in spite of the simple approximations 
and assumptions, the F crit s distribute within a rather narrow portion of total 
F range as shown in Fig. 4. At A = 0.6, the total range of F is < F < 11.8. 
All the F crit distribute within the region 1.5 < F < 4.48, which is roughly 
1/4 of the total region. The region < F crit > ±a occupies only about 11% of 
whole range of F. The similar situation is observed for the other A points. 
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Fig. 4. Phase diagram of the 5-neighbor and 4-state cellular automata. The 
< Pent > /q and the F max jc^ are shown with filled circle. The < F CT n > jc\ + a 
is shown with down triangle and the < F cr i t > /C4 — <r, with up triangle. They are 
joined by lines. A schematic sketch of F m i n in 0.75 < A, and F max in A < 0.25 are 
also shown by the dotted lines. They are drawn based on the arguments in the text. 

The Fig. 4 is still qualitative but rule table space of CA is classified much 
better by using A and F. And it provides us with a deeper understanding of 
the structure of the CA rule table space. The existence of the transition region 
in 0.2 < A < 0.8, means that the edge of the chaos could be found at least in 
this region, which is rather wide range in A. 

Let us discuss the distribution of the rule tables around the transition re- 
gion. All the class IV rule tables are located around transition region as Fig. 
3 shows, but converse is not true. In this region, the three pattern classes 
coexist and the number densities of the class II plus class III rule table are 
considerably larger than that of class IV rule tables. Outside of the transition 
region, the probability of finding the class IV rule table decreases rapidly and 
the pattern classes began to be classified only by the position in the A — F 
plane; in the < F crit > +a < F region, class III rule tables dominates while 
in the F << F crit > —a region, class II or class I CA dominates. The similar 
distributions of the rule tables are observed at A = 0.8, 0.75, 0.7, 0.5, 0.4, 0.3 
and 0.2. 

We proceed to the investigation of the end points of the transition region. In 
the subsection 3.2, we have find the region of A, where both chaotic limit and 
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periodic limit coexist. In the 5-neighbor and K-state CA, it is 1/4 < A < 3/4. 

In A < 0.25 region, not all the group 2 rules could be replaced by the group 
1 rules. Therefore the maximum of N(gl) could not become 256, and it de- 
ceases to zero as A approaches to zero. Then the maximum of F, (F max ) also 
decreases to zero toward A = 0. 

Conversely in A > 0.75 region, not all the group 1 rules could be replaced 
by the group 2 rules. The minimum of N(gl) and therefore the minimum of 
F, (F min ) could not becomes 0. The line F min increases until A = 1. In Fig.4, 
we have schematically shown the F max and F min lines with dotted lines. 

The points where the transition region crosses the F min and F max lines de- 
termine the end points of the coexistence of the two limits, and also the end 
points of existence of the class IV rule tables. We have already found that the 
transition region is found in the region 0.2 < A < 0.8. Therefore the crossing 
points are outside of this region in the case of 5-neighbor 4-state CA. 



5 Conclusions and discussions 



At a fixed A, the patten classes of the CA could be controlled by the numbers 
of the group 1 rules, which has been denoted by N(gl). The maximum of 
N(gl)(= K N ~ l — 1) corresponds to chaotic limit, and N(gl) = , to the 
periodic limit. In the N-neighbor and K-state CA, these two limit exist in the 
region < A < 1 — j^. Therefore in this A region, we could control the patten 
classes by changing N(gl) without fail. The method for it, is the replacements 
of Eq. 8. This provides us with a new method to obtain the rule table of edge 
of chaos. However, we should like to comment that the fluctuation due to the 
random number will make the statements less rigorous. 

This property is studied quantitatively by introducing a new parameter F. 
The maximum of the F corresponds to the chaotic limit and minimum of it, 
the periodic limit. By changing the F using the replacements, we could find 
the region, where the transition of the class II and class III patterns takes 
place. Thus we could obtain the phase diagram of the CA in A — F plane. 

In this article, we have applied the analysis to 5-neighbor and 4-state CA. 
In this case the group 1 rules are further classified into 5 types as shown in 
Table 1, and the phase diagram is obtained in Fig.4. The transition region has 
a rather gentle slope as a function of A, and it extend at least from A = 0.2 
until A = 0.8, which is wider than \ <\<l — \. This explains why the edge 
of chaos is found in the wide range in A for this CA. It may be interesting 
whether or not the slope of the transition region depends on the models. 

In the replacement experiments, we have found the edge of chaos (very long 
transient lengths) in many cases. The examples are shown in Fig. 2. Sometimes 
they are observed in a rather wide range in N 3 or N 2 . This indicates that in 
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many cases, the transitions are second-order like. But the widths in the ranges 
of N 3 or N 2 are different from each other, and there are cases where the widths 
are less than one unit in the replacement of RP2 (first-order like). It is very 
interesting to investigate under what condition the transition becomes first- 
order like or second-order like. The mechanism of the difference in the nature 
of the transition is an open problem and may be studied by taking into account 
effects of group 2, 3, and 4 rules. In these studies another new parameters may 
be found and a more quantitative phase diagram may be obtained. 

It is also interesting to compare the detailed dynamics of class IV CA at 
different points of rule table space in A — F plane by applying Wuencshe's 
method[6] or by using computational mechanics[5]. 

These issues together with finding the points where the transition region 
crosses F max and F min lines in Fig. 4, and the nature of the phase transition 
at these points will be addressed in the forthcoming publications. 
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